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Abstract 



We calculate the long-distance radiative corrections 5 e LD and 5l D to the 
K° — > 7r~e + z/ e and the K° —>■ i\~ [i^v^ decay rates. This analysis includes 
contributions to the long-distance radiative corrections from outside the 
kinematically-allowed three-body Dalitz region and tests the sensitivity of 
the radiative corrections to the hadronic K-tc form factors. A program, 
KLOR, was written to numerically evaluate the radiative corrections and 
-^ ! to generate Monte Carlo events for experimental acceptance studies. The 

K® 3 and the K° 3 long-distance radiative correction parameters are deter- 



mined to be (1.3 ± 0.3)% and (1.9 ± 0.3)%, respectively. We also present 



predictions for the fraction of radiative K® 3 events satisfying various re- 
^ . quirements on final-state photon kinematics. 

O ■ PACS: 13.20.-v, 13.20.Eb, 12.15.Hh 

Oh: 

p; 1 Introduction 

In the Standard Model, the strength and the pattern of weak interactions between 
up- and down-type quarks are encoded in a 3 x 3 unitarity matrix: the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [I]. Years of experimental and theoretical re- 
search have gone into the determination of the elements of this matrix. These mea- 
surements, when combined, test the unitarity of the CKM matrix and thereby deter- 
mine if the Standard Model can fully explain quark mixings. Moreover, measurements 
of the CKM matrix elements can be used to limit proposed extensions of the Standard 
Model. 

The V us element of the CKM matrix is of particular interest because CKM uni- 
tarity is best tested in the first row of the matrix. \V US \ is measured using three-body 
semileptonic decays of the K meson, hyperon decays, and tau decays. Analysis of 
the semileptonic kaon decays, denoted K e3 , leads to the most accurate measure of 
| V us |. Both the charged (ifgj) and the neutral (K^ 3 ) kaon decays provide independent 
measurements of \V US \. 

Up to 2004, analysis of the K^ 3 and the K® 3 decay modes yielded |V^ S | = 0.2200 ± 
0.0026 [2]. The measurements in the K® 3 and Kf 3 modes were in agreement at the 
1% level; however, they were inconsistent with the unitarity of the CKM matrix at 
the 2cr level. This situation soon changed dramatically. A measurement [5] of \V US \ 
from the K^ 3 mode obtained a value of \V US \ that is ~ 3% higher and consistent 
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with CKM unitarity. Results from KTeV [U El E] indicated a value of \V US \, mea- 
sured in the K® 3 and K® 3 modes, that is also consistent with CKM unitarity. The 
precision of the measurement required an understanding of physical effects at better 
than the percent level. Although radiative corrections to these processes had been 
calculated within chiral perturbation theory 00 [9], the KTeV experiment measured 
the angular distributions of radiated photons, which were not modeled satisfactorily 
by currently available treatments [TO]. The present paper was written in response to 
that need. Further measurements of \V US \ were performed by the IHEP |llj . CERN 
NA48 [12], and KLOE [13] Collaborations, leading to a current world average [H] of 
\V US \ = 0.2257 ±0.0021. The KLOE Collaboration utilized radiative corrections doc- 
umented in Ref. [T5]. A concise account of the present status of radiative corrections 
to K® 3 decays is given in Ref. [To] . 

In this paper, we study radiative corrections to the K® 3 and the K® 3 decay modes. 
The total radiative correction to the K® 3 decay rate is described by a parameter 5^. 
This parameter is the sum of two contributions: the long-distance {8 l LD ) and the 
short- distance (5$d) components. The short-distance parameter is known [17J and 
we calculate the long-distance corrections using a phenomenological model originally 
used in Refs. [TBI ESI US EB E21 EH EH- Using the phenomenological model, we 
test the assumption that the variation of the hadronic K-tt form factors is negligible. 
Single-photon radiation from all K® 3 decays are included in our calculation of the 
long-distance radiative correction. A program, named KL0R, was written to numeri- 
cally evaluate the radiative corrections and to generate Monte Carlo events. Monte 
Carlo events may be used by experimentalists to understand detection efficiency For 
example, the detection efficiency of K® 3 decays may be affected by the conversion of 
a radiated photon to an electron-positron pair. 

We also determine the fraction of radiative K® 3 events (normalized by the Kf 3 decay 
rate) that satisfy various constraints on the radiated photon's energy and its proxim- 
ity to the charged lepton. These "radiative fractions" are interesting since they can 
be checked by experiment. Related studies [3 El [9] calculate the radiative corrections 
using chiral perturbation theory (ChPT). 

The paper is organized as follows. In Section [2j we review relevant properties of 
the phenomenological model for the lowest-order contribution to the K® 3 decay. In 
Section El we extend that model to include radiative corrections to the K® 3 decay. In 
Section HI we discuss the numerical method used to calculate the radiative correction. 
In Section El we present results from our calculation of the radiative corrections, while 
Section El concludes. In Appendix |A], we define the relevant loop integrals for this 
analysis and in Appendix [Bl equations used to calculate the virtual matrix element 
are presented. In Appendix E3 we compare our Monte Carlo generator, KL0R, to 
another Monte Carlo generator that can be used to approximate radiative corrections 
in K® 3 decays. 



2 Phenomenological Model of K® 3 Decays 

The K® 3 decay mode and associated kinematics are described by the equation 

K°(p K ) -> 7r-(^) £ + ( Pe ) u t (p v ), (1) 

where £ + (= e + or /x + ) denotes the charged lepton, U£ is the £ neutrino, and px is the 
four-momentum of particle 'X' (= K°, n~ , £ + , or i/^). Conventionally, the lowest- 
order contribution to the K® 3 decay is studied using a phenomenological model |25j . 
Using this model and the momentum parametrization of Eq. (CQ), the lowest-order 
(Born) transition matrix element for the K® 3 decay is [261 I2ZI [28] 

M B = V2G F V us u(p v )^P L v(p e )[f + (t)(p K + P<n)» + (pk - P,)J-(% (2) 

where Gp (= 1.16639 x 1CT 5 GeV) is the Fermi coupling constant, V us is the us 
coupling in the CKM matrix, P L = ~(1 — 7 5 ) is the left-handed projection operator, 
and t = (jpk ~ P-k) 2 - The Mandelstam t variable is the square of the momentum 
transfer to the W + boson (W + — > £ + v). In Eq. (J2J), the functions f±(t) are the K-tt 
hadronic form factors; the form factors arise from the contraction of the quark current 
with the on-shell kaon and pion, 

f+{t)(PK + P«Y + f-(t)(PK-p«r = (ir-( P7r )\sYPLu\K°(p K )). (3) 

2.1 K—tt Hadronic Form Factors 

Before proceeding further, the structure and alternate definitions of the form factors 
are discussed. As shown in Eq. (J2J), the hadronic form factors only depend on t. This 
simple functional dependence is guaranteed by local creation of the lepton pair. Time 
reversal invariance of the K® 3 decay ensures that the form factors are relatively real. 
Instead of using the set of form factors {/+, /_}, an alternate form factor set {/+, 
/o} is typically used. The f+(t) and fo(t) form factors correspond to the vector l - 
and the scalar + exchange amplitudes, respectively. This description of the K-n 
interaction is inspired by a vector dominance model (VDM) [29]. The {/+, /_} form 
factor set is related to the {/+, /o} form factor set by the relation 

/o(t) = /+(*)+ , - -r f-(t), (4) 



m K — m. n 

where rnx and m^ are the K° and n~ masses, respectively. 

Normalized versions of the form factors /+(£) and fo(t) are often analyzed. The 
normalized form factors are defined by the equation, /+,o(£) = /+,o(£)//+,o(0). Re- 
quiring f~(t) to be finite as t — > implies that /+(0) = /o(0); therefore, /+(0) is the 
normalization factor common to both /+(£) and fo(t). We will explicitly write our 
matrix elements in terms of f+ 7 o(t) and /+(0). 

In this paper, we consider three form factor models: the "linear model," the 
"quadratic model," and the "pole model." The linear and the quadratic form factor 



models capture the first- and second-order dependence of the form factor on the 
momentum transfer, t. The normalized linear and quadratic form factor models are 
described by the equations 

/■M)(«) = 1 + A +^> ( 5 ) 

tit 2 
f+,0 (*) = l + A + o — 2 "*" nV ,0~T' 

where the As are model parameters that are measurable by experiment. Since the 
linear and quadratic form factor models are Taylor-series approximations of the "true" 
form factors, they do not describe the higher-order or the large t- dependence of the 
K-tc form factors. 

Unlike the linear and quadratic form factor models, the pole model describes form 
factor ^-dependence beyond second-order and for large momentum transfers. The pole 
model supposes that the behavior of the form factors f+(t) and fo(t) are dominated 
by the exchange of the lightest vector and scalar K* mesons. The pole model expects 
the K* mass poles to be m + ss 890 MeV/c 2 and tuq > m + . In this paper, we consider 
a simple pole model in which the form factors are given by the equation 

/S(«) = Z^-f («) 

m +,0 ~ l 

Modified pole-model schemes (which are not considered in this paper) may include 
the effects of the K* width or the effects of mixing with other vector-meson poles. 

2.2 Born Decay Rate 

Using the matrix element for the K® 3 decay mode [see Eq. (J2J)], the Born contribution 
to the semileptonic decay rate is 



(o) _ ^ F " b K \ v \2 f 2 (f) \ r e m 



where £ (= e or /x), /+(0) is the form factor normalization constant and XL-, is the 
lowest-order phase space integral. The phase space integral is given by the expres- 
sion [30] 

x <») - 4 C dt [ A(t - •*■ m " ) ) ' m t 1 + « ) t 1 - t) 2 ■ (8) 



where rri£ (£ = e or /i) is the mass of the charged lepton, t max = {m 2 K — m 



2 ^2\2 

■k) i 



W-^ UiSrtH^ X®' (9) 

and A (a, b, c) = a 2 + b 2 + c 2 — 2ab — 2bc — 2ca. From Eq. ([HD, one observes that the 
semileptonic decay rate depends on the K-tt form factors. 



3 Radiative Corrections to K® 3 Decays 

The K® 3 decay rate receives corrections from photon emission and exchange. Including 
first-order radiative corrections, the K® 3 decay rate is 

r« 3 = rf £3 + r« 3 = ^IK.| a 4(«) (i + 4) 4). ( 10 ) 

where <5f is the total radiative correction to the Born K® 3 decay rate, /+(0) is the 
form factor normalization constant, and ZjL is the lowest-order phase-space integral 
given in Eq. (jSJ). The radiative correction, 5?, is the sum of two components |31j : 
the short- distance (Ssd) and long-distance (S e LD ) correction. Including the effects of 
quantum chromodynamics (QCD), the short distance radiative correction is JT7] 

Ssd = - 1-- ln-|, (11) 



where a is the fine structure constant, a s = a s (M§) is the strong coupling constant 
at the Z° mass scale, mz is the mass of the Z° boson, and A is the low-energy 
cutoff. The low-energy cutoff is the hadronic mass scale used as a boundary between 
the short- and long-distance loop corrections. The specific cutoff scale is somewhat 
arbitrary and it usually depends on the limitations of the long-distance model; two 
common cutoff scales selected are the proton mass, m p , and the p meson mass, m p . 
The effect of the low-energy cutoff scale is addressed in Sec. EJ 

The long-distance radiative correction is obtained from a model of the interaction 
between the photon, the hadrons, and the charged lepton. Original studies of the 
long-distance radiative correction used an extension of the phenomenological model 
described in Sec. |2] Application of this model to the K® 3 and the Kf 3 decays are in 
Refs. J2DII2T] and Refs. [HI ESI [22] , respectively. More recent studies jZUSE] calculate 
the long-distance contribution to the radiative correction using ChPT. 

3.1 Long Distance Radiative Correction 

We revisit the phenomenological model for radiative corrections to the K® 3 decay 
mode. In this section, we outline the assumptions of the previous studies, discuss our 
modification of these assumptions, and calculate the matrix elements that contribute 
to the K® 3 radiative correction. 

In Refs. [20l EI] , the long-distance radiative corrections to the Kf 3 decay modes 
are computed under the following assumptions: 

AX. Radiative effects are accurately described by first-order perturbation theory in 
a. 

A2. The phenomenological hadronic model [see Sec. [2] characterizes the weak K-it 
vertex. 



4.3. The K-Ti form factors are not modified by the presence of radiative effects. 

4.4. The effect of the pion form factor is negligible. 

.4.5. The t dependence of the K-tt form factors produces a negligible effect on the 
radiative corrections. 

4.6. The experimental apparatus used to detect the radiative K® 3 decays (K® 3 ) 
only detects the charged lepton and the pion; the neutrino and the photon are 
undetected. 

For our calculation, we retain assumptions 4.1-44, we test assumption .45, and we 
remove assumption 4.6 since radiated photons are detected and influence detection 
efficiencies. 

Under the first assumption, 4.1, the Feynman diagrams contributing to the first- 
order radiative corrections are shown in Fig. [U The first-order radiative corrections 
consist of two classes of diagrams: virtual and inner-bremsstrahlung. Virtual dia- 
grams, Fig. [T](a)-(f), involve the emission of a photon by the effective vertex, the 
charged lepton, or the pion and subsequent absorption by the effective vertex, the 
charged lepton, or the pion. The inner-bremsstrahlung diagrams, Fig. [Tfg)-(i), in- 
volve the emission of a real photon by the pion, the charged lepton, or the effective 
vertex. 

Though higher-order corrections in a, electromagnetic corrections to strong- interaction 
graphs, and weak interactions also contribute to the radiative correction, their effects 
are assumed to be small. For studies that model the effects of electromagnetic correc- 
tions to strong-interaction graphs, we refer the reader to the ChPT studies EJ E]. 

To calculate the first-order radiative corrections, we need an expression for the 
"ir°7r~£ + zV effective vertex when at least one of the particle legs is off-shell [see 
Fig. UK a), (g)-(h)]. Under 4.2, the off-shell effective vertex is obtained from the 
on-shell phenomenological model. Therefore, the if°7r - ^ + z^ vertex is 

V K«b, = iy /2G F V us f + (0)P R [(f K +#)/+(0 + (ft -£)/-(0] . (12) 

where p* K and p* are the (potentially) off-shell K° and ir~ four-momenta, t* = (p* K — 
p*) 2 , and the vector j> denotes the contraction of the four-vector p M with the Dirac 
matrices, 7 M . 

In general, one expects Eq. (TT2|) to be an approximation of the true off-shell vertex. 
One potential correction to this vertex is a modification of the K-ii form factors. The 
form factors for the off-shell vertex may depend on kinematical quantities beyond t. 
In fact, using ChPT, Ref. j8] finds that the form factors depend on the additional 
variable s = (p n + Pe) 2 (e.g. f±(t,s)). Since the emission of low-energy photons is 
expected to dominate the radiative corrections, we assume (4.3) that the variation of 
the form factors / +j o is captured by simple ^-dependence. 

The expression for the off-shell K a -K~l + vi vertex may be simplified by using 
energy-momentum conservation and the Dirac equation for the neutrino. The sim- 
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Figure 1: Feynman diagrams for the first-order radiative corrections to the K® 3 decay 
mode. Diagrams (a)-(f) are the virtual corrections while diagrams (g)-(i) are the 
inner-bremsstrahlung corrections. The open circles in diagram (b) and (c) denote the 
self-energy correction to the pion and lepton wavefunctions, respectively. 



plified expression for this vertex is 



V 



K-kIv 



V2G F V us f + (0)P R \2f K f + (t*)-fJ 2 (t*) 



(13) 



where /2(f) = f+(t*) — /-(£*) is the difference between the /+ and /_ form factors. 

We assume that the electromagnetic interaction with the pion is point-like (.4.4). 
Deviations from the point-like assumption are captured by the introduction of pion 
form factors. For experimental studies of the strength of pion form factors consult 
Ref. [H]; in this study, we assume pion form factors are negligible for current exper- 
imental sensitivities. Previous studies (22] of radiative corrections to the K^ 3 decay 
mode include the effect of the pion form factors. Though the long-distance radiative 
correction to the K& mode from Ref. [22] is larger in magnitude (5ld = —2%) than 
other studies that assume a point-like pion (Sld = —0.45%) [HUE], the discrepancy 
is believed to be calculational rather than physical. 

In Sec. 13.1.11 and Sec. 13.1.21 the matrix elements for the inner-bremsstrahlung 
contribution and the virtual contribution to the radiative corrections are calculated. 
The numerical methods used to evaluate the inner-bremsstrahlung and virtual matrix 
elements are discussed in Sec. SJ 



3.1.1 Inner-Bremsstrahlung Contribution 

Previous studies of single-photon radiation in K® 3 decays using the aforementioned 
phenomenological model are found in Refs. J2QI I2D 1231 123- In Refs - Effl 121], the 
K-7T form factors were taken to be constant (.4.5) and the final state neutrino-photon 
system was assumed to be "undetected" (.46). Since the sensitivities of modern 
experiments are expected to probe the form factor dependence of K® 3 decays, we 
relax .45. Current experiments also measure the radiated photon from K® 3 decays; 
therefore, we remove .46 in our analysis. 

In Refs. [23l [23], the K® 3 matrix element was derived assuming a linear form 
factor model and using soft-photon theorems [321 S3]- Unlike Refs. [201 EI], their 
study does not assume that the photon is undetected; therefore, radiative photon 
spectra may be generated for use in detector-acceptance studies. In this section, we 
derive the inner-bremsstrahlung matrix element for the quadratic and the pole form 
factor models and we revisit the constant and linear models. 

The inner-bremsstrahlung correction to the K® 3 decay rate consists of single pho- 
ton radiation from the pion, the lepton, and the effective vertex [see Fig. [Tfg)-(i)]. 
The kinematics for the inner-bremsstrahlung diagrams are described by the decay 
equation, 

K°( PK ) - 7r-(p w ) £ + ( P t) v t {p v ) 7 (*), (14) 

where k is the photon four-momentum. 

Using the off-shell effective vertex, Eq. (fl3l . the contribution to the inner-brems- 
strahlung matrix element from radiation off the pion (Fig. 11(g)) and the lepton 
(Fig.[H(h))is 



M (9)+W = eV2G F V u J + {Q,)u{p v )P R \\2i> K U{t 1 )+mJ 2 {t l 
tyicf+fe) + mtf 2 (b 



e*(k) -p % 



e*(k)-p i + ^(k) 



k ■ Vf 2k ■ 



k-Pn 

+ Mt2)f(k)\v(pt), (15) 



where e^(k) is the photon polarization vector, me is the mass of the charged lepton, 
and f2(t) = /+(£) — f-(t). In Eq. (fl5l) . the form factors are evaluated for two different 
values of the momentum transfer to the W + gauge boson: t\ and t 2 - When the photon 
is radiated from the pion, the momentum transfer is t\ = (pk — Pit — k) 2 = (pe + p u ) 2 
and when the photon is radiated from the charged lepton the momentum transfer is 

h = (Pk ~ Pn) 2 - 

The last component of the inner-bremsstrahlung matrix element is from photon 
radiation off the effective vertex [see Fig. QJi)]. The expression for the u K Q ir~l + v^ n 
vertex is obtained by requiring the matrix element to be gauge invariant. In Ref. [2D1 
W\\ the K°iT~e, + vi r ) vertex is obtained by applying the principle of "minimal coupling" 
to the weak-interaction Lagrangian, while Refs. [231 123] uses the Ward identity [33] 
to determine the K°TT~£ + i , g'-f vertex. We illustrate both of these methods and we 
use the latter to determine the forms of the K°t^~1 + v^ vertex for non-constant form 
factors. 
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Assuming the form factors are constant, the weak- interaction Lagrangian is 

C w = V2G F V US fa-fP^i [(/+ + f-Hird^K + (f + - f_)<p K dM , (16) 

where tp and <fi denote the fermion and complex scalar fields, respectively. Replacing 
the partial derivative with the covariant derivative (minimal coupling), 

d^ .-► D^ = d„ - ieQAp, (17) 

where A^ is the photon field, one obtains the K°7r~£ + i>£'-f effective vertex, 



V™" 7 '" = -iV8^cG F V us P R f 2 r. (18) 

Combining all three inner-bremsstrahlung graphs, the matrix element for the constant 
form factors is 



M h JZ = eV2G F V us f + (0) u{p u )P R [2j> K + 



mt\ 



e*(k)-p n e*(k)- Pl jjg(fe) 
k ■ Pit k-pi 2k- Pi 



v{pi 



^ (19) 
Since f +y2 are normalized and independent of t they have been set to 1. Taking 
e M (/c) — > k^ in Eq. (TT91 . one observes that M^ns ~^ an d therefore the Ward 
identity is satisfied. 

Extraction of the fT°7r~£ + z^7 vertex through minimal coupling works well for the 
constant form factor model. When form factors are no longer constant and when the 
matrix element depends on two different values of t (e.g. t\ and t 2 ), this method is 
insufficient. For this scenario, we determine the K°ir~£ + i>£'y vertex by requiring the 
inner-bremsstrahlung matrix element to be zero when the photon polarization vector 
e tl (k) is replaced by the photon four- momentum (Ward identity). This method does 
not uniquely determine the K°ii~£ + vi r ) vertex; the vertex expression may differ by 
structure- dependent terms. For a discussion of structure-dependent terms we refer 
the reader to Ref. [231 123]- 

Using this methodology, the i^°7r~£ + z/^7 vertex for the linear form factor model is 

v K«*n, „ = _ 2v /g^ G F V us f + (0)P R {f 2 {h) Y ~ 2/;(0) [2j> K - U (p K - p,)"} , 

(20) 
where /+(0) (= A + /m^) is the slope of the normalized form factor. The vertex 
interaction for the quadratic model is 

VfS 1 '" = -*v/8^ G F V us f + (0)P R {f 2 (t 2 ) ^ 

-2 [2} K -h\ [/;(o) + |(t! + t 2 )K(oj\ (pk-p.Y (2i) 

- Mm 2 K - ml) [f»(0) - f»(0)} (p K - p«r} , 
where f' + (0) (= X' + Q /m%) and f" (0) (= A" o/ m ^)- The vertex interaction for the 



pole model is 



v; 



p^if 7 ' " = -*VS^ G F V us f + (0)P R {f 2 (t 2 ) 

-2 \2i> K - fa] D + {t u t 2 ){ PK - Pn r 



r 



-2fa(m 2 K - ml) 



m i — rrir 



'+"'0 



M (p* - ?*)" 



(22) 



x[/o(ii)^(ti,t 2 ) + / + (t2)A(ti,t 2 )]} 
where m + and mo are the poles of the K* mesons and 

D +>0 (t 1 ,t 2 ) = m; 2 / +i0 (ti)/+,o(t 2 )- 



(23) 



Combining the K°tt i^v^ contribution with Eq. (Tl5|) . the inner-bremstrahlung 
matrix element is 



M hiem = eV2G i ,y us / + (0)«(^)P R U 2 ^/+(*i)H-m,/ 2 (t 1 
- 'ty K f+(t 2 )+mJ 2 (t 2 



e*{k)- Vl+ mk) 



k 



2k- 



e*{k) -p n 
k-v-n 

+ X}v( Pi ), 



(24) 



where X is non-zero for non-constant form factor models. The three following equa- 
tions define X for the linear, the quadratic, and the pole model: 



X 



lin 



2£/'(0)[2 2 J Jf -& 



* quad = 2£ [2jfe - fa] [/;(0) + |(ti + t 2 )/:(0)] - £&(m£ - "4) [f'm - /o'(0) 
A'poie = 2£[2i> K - fa] D+it^h) 

+2&(™ 2 k ~ ml) i 11 ^] [fo(ti)D + (t h t 2 ) + U{t 2 )D {t u t 2 )} , 



(25) 
where £ = e*(fc) ■ (j>k — Pn), and m + and m are the vector and scalar K* masses of 
the pole form factor model. 

The terms in Eq. ff25l) are necessary to maintain the gauge invariance of the matrix 
element. In our numerical evaluation of the inner-bremsstrahlung matrix element [see 
Sec. Sj, we introduce a small photon mass A to regulate the infrared divergence. If 
the contribution from Fig. [T](i) to the matrix element was not included, then the 
matrix element would grow like ~ E 1 /\. This behavior arises from the longitudinal 
component of the photon polarization vector in the term f 2 (t 2 )j£(k) and also from the 
the fact that ti ^ t 2 . Including the effect of Fig. [T](i) in the matrix element, the Ward 
identity is maintained and the matrix element has the correct behavior as the photon 
mass, A, goes to zero. 
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3.1.2 Virtual Contribution 

The first-order virtual corrections consist of three types of diagrams: self-energy 
corrections (Fig. IU(b)-(c)), inter-particle photon exchange (Fig. [Qa)), and effective 
vertex interactions (Fig. [Tfd)-(i)). Previous studies [201 E] calculate the virtual 
corrections assuming that the form factors are constant. In this section, we recalculate 
the virtual matrix element assuming constant form factors and we investigate two 
methods to calculate the virtual contribution for non-constant form factors. 

Self-energy corrections to the pion and to the charged lepton are obtained by wave- 
function renormalization of the weak Lagrangian [see Eq. (TToT) ]. The wavefunction 
renormalizations for the pion and the charged lepton are 



SZ e 

bZ„ 



An 



2Bi(m|; mj, A 2 ) + Am\B' x {m\\ mj, A 2 ) + 8m 2 e B' {mj; mj, A 2 )J (26) 
1 + 2B (ml ; m 2 , A 2 ) + (4m 2 - A 2 )^(m 2 ; m 2 , A 2 )l , (27) 



where B(-) and -B'(') are the two-point loop integral [12] and its derivative, 
B'(p\, m 2 ,, m\) = dB(p\, m§, m 2 )/<9p 2 . The loop integrals are calculated using a cutoff 
regularization and are defined in Appendix [A] The self-energy contribution to the 
virtual matrix element is 

M {b)+{c) = ^{SZe + SZ^M 3 , (28) 

where M. B is the lowest-order K® 3 matrix element [see Eq. (J2J)]. 

The inter-particle photon exchange diagram is calculated using Eq. (TI3"1) for the 
K°7t~£ + b>£ vertex. Since the form factors are assumed constant, they may be pulled 
outside of the photon four-momentum integration. For non-constant form factor 
models, the form factors will depend on the photon four-momentum and directly 
modify the photonic loop integration. 

For the constant form factor model, the vertex correction consists of the dia- 
grams (d) and (e) of Fig. [TJ Since the minimal coupling anzatz does not require a 
i^°7r~£ + z/£77 vertex, Fig. \Hf) is neglected in the virtual correction calculation. The 
i^°7r~£ + z/£7 vertex from Eq. f[T^]) is used to calculate the effective vertex interaction 
graphs. As with the evaluation of the inter-particle photon exchange diagram, the 
constant form factor does not modify the integral over the photon four-momenta. 

Combining the contributions from each of the virtual diagrams, the virtual matrix 
element for a constant form factor model is 

A^irt = a ^ GfVu j + ( 0) u {p v )P R \Ai> K - Bfa } v(p £ ), (29) 

47T 

where A and B may be written in terms of photon loop integrals. A and B are given 
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by the equations 



A = f + [1 + 4B (m£; ml A 2 ) + 25! K 2 ; mj, A 2 ) 
+ (4m 2 -A 2 ) J B (m 2 ;m 2 ,A 2 ) 
+8mjB' (mj; mj, A 2 ) + AmjB[(mj] mj, A 2 ) 
-8x rf C + 4m 2 d - 4 (2x rf + m 2 ) C 2 ] 
+/a [2Bo(sir/; ""#, mf ) - 5 (m 2 ; m 2 , A 2 ) 
+4B (mj; mj, A 2 ) + 2B l (mj; mf; A 2 ) 
+Bi (m 2 ; m 2 , A 2 ) + S x (m 2 ; A 2 , m 2 ) 

»/ 2 C 2 ] 



(30) 



+2A 2 C -4x^C 1 + 2m 2 C 2 



and 



B= f + [-4s ni C 2 ] 

+/ 2 [| + 2E (^; m 2 , m 2 ) + E (m 2 ; A 2 , m 2 ) 
-45 (m 2 ; mj, A 2 ) + ^(m 2 ; m 2 , A 2 ) 
+B 1 (m 2 ; A 2 , m 2 ) - B 1 {mj; mj, A 2 ) 

+4m 2 5 (m 2 ;m 2 ,A 2 ) (31) 

+ ±(4m 2 -A 2 ) J B (m 2 ;m 2 ,A 2 ) 
+2mjB[(mj;mj,X 2 ) 
+ (2A 2 - 4x rf )C + (2m 2 - Ax^)C x 
+ (2m 2 - 4x n£ )C 2 } 

where s^- = (p, + p,,) 2 , a?y = p* • p,-, and C, = C^m 2 , s^, m 2 ; A 2 , m 2 , m 2 ). 

In Eq. (I30 l -( 13T1) . a small photon mass A is introduced to regulate the infrared (IR) 
divergence. Since the K°ir~£ + Uf vertex is an effective interaction, the interaction is 
not renormalizable and will depend on an ultraviolet (UV) cutoff scale, A. Assuming 
a constant form factor model, the virtual matrix element is logarithmically divergent 
(.M virt ~ In A 2 ). A cutoff regulator is used when calculating the loop integrals; we 
select the UV cutoff such that A > m n , m,£. We take A to be the proton mass (m p ). 
In Sec. El we discuss the sensitivity of the radiative correction to the UV cutoff scale. 

For non-constant form factors, Eq. (I30~1) - (I3~TI) will be modified. In this paper, we 
use two methods to study the effects of non-constant form factors in the virtual calcu- 
lation. The first method, Method I, is an approximate method in which the constant 
form factors / +> o in Eq. (I30l-(l311 are replaced by f +j0 (t), where t = (px — Pn) 2 - 
Using this method, each K® 3 event is modified based on the size of the momentum 
transfer to the lepton-neutrino system. 
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The second method, Method II, assumes the pole form factor model and it in- 
cludes the effect of the photon four-momenta on the form factor. For instance, in 
the inter-particle photon exchange diagram (Fig. [0(a)), the form factor depends on 
the momentum transfer t\ = (px — P-k — k) 2 . Therefore, the /c-dependence of the 
form factor will modify the photonic loop integration. Expressions for A and B in 
the Method II are in Appendix [BJ For Method II, we also use cutoff regularization; 
therefore, a new UV cutoff scale is selected that is larger than the mass of the K* 
poles. 

Results from both methods are presented in Sec. [5j 

4 Numerical Methods 

To determine a numerical value for the long-distance radiative correction, 5 e LD , we 
integrate the squared, summed inner-bremsstrahlung and virtual matrix elements 
over their respective final-state phase spaces. To aid in this evaluation, we developed 
a program named KLOR (Kaon Leading Order Radiation). Using KLOR, we numerically 
evaluate the inner-bremsstrahlung and the virtual matrix elements for both constant 
and non-constant form factor models. In this section, we discuss the methodology of 
the numerical calculation. 

Using the expressions for the inner-bremsstrahlung and the virtual matrix ele- 
ments in Sec. [3] KLOR numerically squares each matrix element and sums each of 
them over polarization- and spin- states. The squared, summed matrix elements are 
integrated over phase space by two methods: Vegas integration [35] and traditional 
Monte Carlo integration [36]. The Vegas integration is used to evaluate the total 
decay rate, while the traditional Monte Carlo integration is used to check the Ve- 
gas integration and to evaluate the decay rate when final-state requirements on the 
photon kinematics (e.g. angular and/or energy constraints) are imposed. An op- 
timized, phase-space sampling technique is used to efficiently generate unweighted 
event samples. 

Since our calculation of the decay rate is numerical, proper regulation of the IR 
divergence is necessary to obtain finite results. The IR divergence is regulated by 
introducing a small photon mass A. For the inner-bremsstrahlung calculation, the 
photon mass limits the kinematically accessible four-body phase space. Since the 
photon is massive, one needs to sum over the longitudinal polarization in addition to 
the two transverse polarizations. In the virtual calculation, the photon mass modifies 
the loop integration over the photon momentum. Tests of the numerical stability 
of our IR regularization scheme show that KLOR's results are numerically stable for 
photon masses from ~ 10 -5 to 10~ 14 GeV/c 2 . 

The phenomenological K TT~£ + ug interaction is effective (non-renormalizable); there- 
fore, this model describes the interaction up to a UV cutoff scale, A. We take the UV 
cutoff, A, to be the low-energy cutoff described in Sec. [3J recall that the low-energy 
cutoff is the boundary between the long- and short- distance corrections. The UV 
cutoff is used to regulate the loop integration - cutoff regularization. In particular, 
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the photon energy is integrated up to the UV scale A. 

Loop integrals in the virtual matrix element are computed using the numerical 
package Looptools J37J EH]- Looptools calculates the value of loop integrals using 
analytic expressions derived for dimensional regularization. As long as the cutoff scale 
is taken to be large; evaluation of the loop integrals in the cutoff regularization can 
be related to the corresponding expressions in dimensional regularization with minor 
modifications. 

In Appendix [Cj we compare results from KLOR to other approximate techniques. 
In particular, we analyze distributions produced by KLOR and another program called 
PHOTOS [10]. PHOTOS is a Monte Carlo event generator that can be used to approximate 
the radiative corrections to the K® 3 decay modes. 

5 Results 

In Sec. [31 we developed a phenomenological framework to describe radiative correc- 
tions to the K® 3 decay modes and in Sec. HI we described a program, KLOR, that 
was written to implement this framework. In this section, we use KLOR to illustrate 
the structure of the corrections, to determine the 5 e LD parameter, and to calculate 
experimentally measurable observables of the model. 

Before delving into the aforementioned topics, we would like to illustrate the size 
of the radiative corrections. One approach is to examine the invariant mass of the 
pion/lepton system, denoted m % £. In Fig. [21 we plot the Born and the next-to- 
leading order (NLO) differential decay rates as a function of m ne (Fig. 121(a)) and m nfl 
(Fig. 12(b)). In the notation of Eq. (ITOl . the Born decay rate is T ( K ' i3 and the NLO 
decay rate is T { K ' t3 + T { K ' e3 . 

The effect of radiative corrections on the Born decay rate is illustrated by the shifts 
of the NLO distribution from the Born distribution. In Fig. [2(c) and Fig. [2(d), we 
plot the ratio of the NLO and the Born m n e distributions for the K® 3 and K® 3 modes, 
respectively. In each of these plots, we assume a pole model [see Eq. ([6])] for the 
hadronic K-tc form factors. The parameters of the pole model, m + = 0.882 GeV/c 2 
and m = 1.174 GeV/c 2 are from a recent measurement by the KTeV Collabora- 
tion [B|. 

As observed in Fig. [21 radiative corrections do not uniformly modify the m ne 
distribution. Radiative corrections enhance the Born decay rate for m^ masses below 
0.3 GeV/c 2 and suppress the Born decay rate for m ne masses above 0.3 GeV/c 2 . A 
typical correction to the m^ distribution is ~ 5% but can be greater for m^ e masses 
near the kinematic boundaries. Radiative corrections to the m nfl distribution tend 
to be more uniform and generally enhance the decay rate. A typical K® 3 correction 
to the m^ distribution is ~ 3% but can be larger near the endpoints of the m nfl 
distribution. 

Now that we have a basic understanding of the size of the radiative corrections in 
the K® 3 and the K® 3 modes, we investigate the structure of the radiative corrections. 
In particular, we study radiative corrections over the K® 3 and the K Q 3 three-body 

14 




hi e (GeV/c") 



m^ (GeV/c z ) 



Figure 2: (a) The invariant mass distribution of the pion-electron system (m e7r ) using 
the Born (solid line) and the NLO (dashed line) matrix elements, (b) The invariant 
mass distribution of the pion-muon system (m M7r ) using the Born (solid line) and the 
NLO (dashed line) matrix elements, (c) and (d) are the ratio of the NLO and the 
Born distributions for the m en and m en distributions, respectively. 

Dalitz regions. 



5.1 Dalitz Distribution of the Radiative Corrections 



In Fig. [3](a) and Fig. IH(a), we plot the lowest-order decay rates for the K® 3 and 
the K° 3 decay modes over their respective three-body Dalitz regions. These plots 
illustrate the dependence of the Born decay rate on the reduced energies 

2-C/7T 1Ep 

y = — , z = — e - } 32 

rax rax 

where E n (Eg) is the pion (charged lepton) energy in the kaon center of mass. 

When radiative effects are included, the Born Dalitz distributions are modified. 
Using Method I [see Sec. I3.1.2J . the NLO decay rate for the K® 3 and K® 3 modes are 
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Figure 3: (a) Dalitz plot of the Born decay rate for the K® 3 mode, (b) Dalitz plot for 
the K® 3 decay rate including radiative corrections. Plot (b) uses the pole form factor 
model (m+ = 0.882 GeV/c 2 , m = 1.174 GeV/c 2 ) [6] and Method I to calculate the 
virtual diagrams. The cutoff scale A = m p is used. 



plotted in Fig. EJh) and Fig. IU(b), respectively. In general, radiative corrections mod- 
ify the Born decay rate both inside and outside the Dalitz region. Inside the Dalitz 
region, both virtual and inner-bremsstrahlung diagrams contribute to the NLO decay 
rate; while outside the Dalitz region, only inner-bremsstrahlung diagrams contribute. 
The inner-bremsstrahlung contribution to the K® 3 decay rate outside the Dalitz re- 
gion is more pronounced in the K® 3 decay mode than in the K° 3 decay mode since 
there is a higher probability of radiation from the electron. 

The long-distance radiative correction to each point in the three-body Dalitz plane, 
denoted 8 l LD {y : z) : may be obtained by taking the ratio of the NLO and the Born 
distributions and then subtracting unity. S e LD (y, z) is interesting since it describes the 
structure of the radiative correction; for a given (E^, Eg), 8 e LD (y, z) characterizes how 
the Born distribution is modified. 

In Fig. El we plot S LD (y, z) using Method I for both the K® 3 and the K® 3 modes. 
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Figure 4: (a) Dalitz plot of the Born decay rate for the K® 3 mode, (b) Dalitz plot for 
the K Q 3 decay rate including the radiative corrections. Plot (b) uses the pole form 
factor model (m+ = 0.882 GeV/c 2 , m = 1.174 GeV/c 2 ) [6] and Method I to calculate 
the virtual diagrams. The cutoff scale A = m p is used. 



Radiative corrections to the K® 3 Dalitz region can be significant, ranging from —15% 



to 



15%. Radiative corrections to the K® 3 



Dalitz region are smaller than for the 
K® 3 mode; they range from —2% to 8%. Since there is no contribution to the Born 
decay rate from outside the Dalitz region, S e LD (y, z) is infinite in this region. In Fig. [5j 
we have set 5 e LD (y,z) outside the Dalitz region to zero. Depending on the technique 
used to measure the K® 3 decay rate, the contribution from outside the Dalitz region 
may need to be included in the S LD calculation. In Sec. 15. 2\ we address this and other 
issues affecting the 5 l LD calculation. 

Before we calculate 8 LD , we compare the two methods, Method I and Method II, 
proposed in Sec. 13.1.21 Recall that Method I and Method II are two ways to introduce 
non-constant form factors into the virtual matrix element calculation. We compare 
these methods by generating radiative distributions and examining the differences 
between them. 
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Figure 5: Long-distance radiative correction, 5 e LD , to K® 3 (a) and the K® 3 (b) Dalitz 



plots. Both plots use the pole form factor model 



777m 



0.882 MeV/c 2 , m = 1.174 



MeV/c 2 ) [6] and Method I to calculate the virtual diagrams. The UV cutoff (A) for 
the long-distance calculation is taken to be the proton mass. 



In Fig. [6](a), we plot the difference between the NLO and Born decay rate as a 
function ofm^,,. The solid (dashed) line is produced using Method I (Method II). The 
m m distribution using Method II is shifted higher than that produced by Method I; 
however, the size and structure of the distributions are similar. 

In Fig. EKb), we plot the y (= 2£' e /m ft -) distribution for the difference between 
the NLO and Born decay rates in the 0.7 < z < 0.8 "band" of the Dalitz plane. As 
with the m ne distribution, the radiative correction from the Method I and Method II 
have very similar sizes and shapes. Inspection of other bands in the Dalitz plane lead 
to similar findings. 

Since the radiative correction distributions produced by Method I and Method II 
produce commensurate spectra, one expects that the distribution of K® 3 events pro- 
duced by Method I and Method II to be identical for current experimental precision. 
The distributions from Method II are in general more positive than those of Method 
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Figure 6: (a) The m^ distribution of the difference between the NLO and the Born 
K® 3 decay rates. Method I (solid) and Method II (dashed) are presented, (b) The y 
distribution of the difference between the NLO and the Born K® 3 decay rates in the 
0.7 < z < 0.8 band of the K% Dalitz plane. Method I (solid) and Method II (dashed) 
are presented. 

I; therefore, we expect that the total radiative correction for Method II to be slightly 
larger. 



5.2 



The 8 LD Parameter 



The total radiative correction, denoted <5f , to the Born K® 3 decay rate consists of a 
short- and long-distance component. Using the phenomenological model [see Sec. [3], 
we calculate 6 e LD and estimate its uncertainty. Before we can determine S e LD , we must 
first specify the analysis technique used by experimentalists to measure the K% decay 
rate. 

In previous studies [201 ED E] , the authors assume that the experimental apparatus 
only measures the charged lepton and the pion; the neutrino and (if present) the 
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photon are "invisible." The observed pion and charged lepton momenta are then fit to 
three-body kinematics with zero missing mass. If a Kf 3 event produces a '#7r' signature 
that lies outside the three-body Dalitz region, then it is excluded from the decay rate 
measurement. Note that this scenario only works for an ideal experiment for which 
the kaon energy is known and where the analysis cuts are such that events can not 
migrate from outside the three-body Dalitz region. For this experimental analysis 
technique, the long-distance radiative correction parameter is the total correction to 
the Born decay rate inside the three-body K® 3 Dalitz region. In this scenario, we 
denote the long-distance radiative correction 5 e LD . 

We consider an experimentally motivated approach to the analysis in which the 
experiment is sensitive to and detects radiated photons in the final state. There- 
fore, an analysis of the K® 3 decay rate will include radiated photons. Moreover, we 
also assume the experimental analysis includes K® 3 signatures with measured pion 
and charged lepton energies outside the three-body Dalitz region. In this analysis, 
5 e LD includes contributions from both inside and outside the three-body Dalitz plot 
[see Fig. [3]. Note that since we perform a numerical analysis [see Sec. H], one can 
compute S e LD , which excludes events from outside the three-body Dalitz region. 

In addition to the experimental technique, 8 LD will depend on the cutoff scale A. 
As described in Sec. [3j both the short- and long-distance radiative corrections depend 
on the scale A which denotes the "boundary" between the short- and long-distance 
loop corrections. Below A, radiative corrections to the K® 3 decay rate are described by 
the phenomenological model; while, above A, radiative corrections are described by 
the Standard Model. In our analysis, we associate uncertainty in the total radiative 
correction (Sj>) due to the cutoff scale choice with uncertainty in the long-distance 
radiative correction. 

The choice of the cutoff scale is somewhat arbitrary and typically depends on the 
"physics" and the calculational limitations of the long-distance model. Since we use 
the cutoff regulator to compute the loop integrals, we take A to be larger than the 
masses appearing in the loop integrals. In Method I and Method II, we take A to be 
m p and 3 GeV/c 2 , respectively. Though the assumption that the long-distance model 
describes physics at the cutoff scale A = 3 GeV/c 2 is likely optimistic, integration of 
the loop integrals, which depend on the K* poles, requires a higher cutoff scale. 

In Table [TJ and Table [2J we compute the short- and long-distance radiative cor- 
rections for various cutoff scales ranging from m p (= 0.7711 GeV/c 2 ) to 3 GeV/c 2 . 
For Method I, changing the cutoff scale from m p to 3 GeV/c 2 results in changes to 
5 e LD , 5l D , and to 5sd of ~ 0.7%; however, 6? (= Ssd + ^ld) remains stable to 
within ~ 0.2%. Therefore, a reasonable uncertainty in 5 e LD due to the cutoff scale 
uncertainty is 0.2%. A similar inspection of the K° 3 mode indicates an uncertainty 
to 81d fr° m the cutoff scale is 0.1%. 

In addition to the uncertainty associated with the cutoff scale ambiguity, there is 
an uncertainty in the radiative corrections due to the method used to calculate the 
virtual corrections. Using Method II, the total radiative corrections to the K® 3 and 
the K° 3 modes are slightly larger than the corresponding values obtained via Method 
I. Therefore, for 8 e LD and 5^ D we take an additional uncertainty of 0.1% and 0.2% due 
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Table 1: Long- and short-distance [3QJ radiative corrections to the K® 3 decay modes. 
Results are presented for different virtual correction methods (I and II), for different 
form factor models, and for various UV cutoff scales. We use the measured values of 
the form factor parameters taken from Ref. [Bj. 



FF Model 


A 


5sd 


Sin 


5 e T 




(GeV/c 2 ) 


(%) 


(%) 


(%) 


Method I 










Constant 


m p 


2.27 


1.26 


3.53 


Constant 


m p 


2.17 


1.31 


3.48 


Constant 


2m p 


1.94 


1.55 


3.49 


Constant 


3 


1.64 


1.73 


3.37 


Linear 


m p 


2.17 


1.32 


3.49 


Quadratic 


m p 


2.17 


1.33 


3.50 


Pole 


m p 


2.17 


1.33 


3.50 


Method II 










Pole 


m p 


2.17 


1.43 


3.60 


Pole 


1m p 


1.94 


1.84 


3.78 


Pole 


3 


1.64 


2.11 


3.75 



to the method that is used to calculate the virtual corrections. 

Therefore, after inspection of the shifts to the total radiative corrections in Tabled] 
and Table [2j we take the uncertainty on the long-distance radiative corrections due 
to the cutoff scale and the virtual calculation method to be 0.3% for both 5 e LD and 

°LD- ___ 

Using Method I and Method II [see Sec. 13.1.2] . we also study the sensitivity of 
the radiative corrections to variations in the form factor parameters. In Fig. [7J, we 
use Method II to calculate the total radiative corrections 5^ and <5j> as a function 
of the pole model parameter m+; the pole model parameter mo is held constant in 
these plots. In Fig. [3(a), the pole model parameter m = 1.174 and in Fig. [7(b), 
we generate the curves using m = 1.174 GeV/c 2 and m = 0.9 GeV/c 2 . For the 
K® 3 decay mode, 5 e LD is insensitive to m since the f (t) form factor is multiplied by 
the square of the electron mass. 

Using Method I, we find that the total radiative corrections to the K® 3 and 
K® 3 modes are insensitive to changes in the m+. For Method I, the insensitivity 
of 5 e LD and 51 D to changes in the form factor parameters is likely a result of the 
simplistic method used to introduce t dependence of the form factors into the virtual 
correction, / +0 —> f+,o(t) [see Sec. 13.1. 2j . 

Using Method II, the total radiative correction shows variation with changes in 
the pole model parameters. Using this method the total radiative correction varies 
~ 0.3% for the K® 3 and ~ 0.2% for the K ^ modes over the m + range 0.6-1.5 GeV/c 2 . 
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Table 2: Long- and short- distance [3H] radiative corrections to the K® 3 decay modes. 
Results are presented for different virtual correction methods (I and II), for different 
form factor models, and for various UV cutoff scales. We use the measured values of 
the form factor parameters taken from Ref. [Bj. 



FF Model 


A 


5sd 


Kn 


6£ 




(GeV/c 2 ) 


(%) 


(%) 


(%) 


Method I 










Constant 


m p 


2.27 


1.79 


4.06 


Constant 


m p 


2.17 


1.88 


4.05 


Constant 


2m p 


1.94 


2.17 


4.11 


Constant 


3 


1.64 


2.35 


3.99 


Linear 


m p 


2.17 


1.93 


4.10 


Quadratic 


m p 


2.17 


1.93 


4.10 


Pole 


m p 


2.17 


1.92 


4.09 


Method II 










Pole 


m p 


2.17 


1.96 


4.13 


Pole 


1m p 


1.94 


2.38 


4.32 


Pole 


3 


1.64 


2.66 


4.30 



Though Method II indicates variation of 5 e LD and Si D with changes in the form factor 
parameters, for reasonable changes in the form factor parameters the variation to the 
radiative correction parameters are small (~ 0.04%). Therefore, the assumption that 
the variation of the K-tc form factor parameters produce a negligible effect on the 
radiative corrections is reasonable. 

Using the phenomenological model and a cutoff scale of A = m p , we find that 
the long-distance radiative corrections to the K® 3 and the K® 3 decay modes are 1.3% 
and 1.9%, respectively. For this model, the uncertainties on 5 e LD and <5^ D result from 
the cutoff scale/virtual calculation method (0.3%) and the form factor parameter 
dependence (0.04%). Therefore, we take the uncertainty on 5 e LD and 5^ D to be ~ 0.3%. 

In previous studies of the long-distance radiative corrections, Refs. [201 EI] found 
that 5 e LD = 1.5% and 51 D = 2.0% [40]. Recall that 5 e LD (£ = e, or /x) corresponds to 
the long-distance radiative correction inside the K® 3 three-body Dalitz regions. Since 
$ld < S e LD , one infers that 5 e LD > 1.5% and 5£ D > 2.0% from Ref. [201 [21] • Using our 
calculation with a constant form factor model, 5 e LD = 1.3% and 8^ D = 1.9% [see Ta- 
bled] and Table [5]. The cause of the discrepancy between the S e LD values is unknown; 
however, Ref. [S] has noted a formerly undocumented error in the K^ 3 calculation 
of Ref. [201 [21] • Therefore, one possible explanation is an unknown error in the 



K% calculation of Ref. [201 EH- 



A more recent study [8] of the K® 3 decay rate using ChPT finds that for A 



rn 



p- 



5 e LD is 0.6%. When the cutoff scale is taken to be m p , 5 e LD should increase by ~ 0.1%. 
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Figure 7: (a 

function of the K* pole mass m + . Curve is generated using Method II and ttiq = 1.174 
GeV/c 2 . (b) Total radiative correction, 8? (= S^ D + Ssd), to the K® 3 decay rate as 
a function of the K* pole mass m + . Curves are generated using Method II and 
m = 1.174 GeV/c 2 (solid) and m = 0.9 GeV/c 2 (dashed). 



In our calculation, we find that for the K® 3 decay rate, the contribution to 5 e LD from 
outside the K® 3 Dalitz region is ~ +0.5%. Therefore, we infer that 8 e LD ~ 1.2% from 
Ref. [8J, which agrees with our calculation. 

Predictions of the phenomenological model may be tested experimentally. One 
way to test the long-distance radiative corrections 5 e LD and S^ D is to test the ratio 



T t 



1 + 8^ 



LD 



LD 



1 + 8' 



(33) 



LD 



Using our radiative correction parameters, we find that Tld = 1-0058 ± 0.001. The 
uncertainty in J~ld is obtained from Table Q] and Table EJ note that the common 
cutoff scale uncertainties likely cancel in the ratio. Using Eq. (TTOl) and assuming 
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lepton universality, the ratio Tld equals 

(34) 



r ^ 3 jfofr) 



\:.-i 



W 



where Y K o /T K o is the ratio of the K® 3 and the K® 3 decay rates, I? J\i) / I(L(\i) is 
the ratio of the lowest order K® 3 and K® 3 phase space integrals [see Eq. (jSJ)], and the 
Aj are form factor model parameters. 

Using the measured decay rates and the form factors from Ref . |I] , we find that the 
experimental value of Tld is 1.0027 ± 0.0048. This is consistent with our calculated 
value. 

5.3 Radiative Fractions 

In addition to the fraction Tld defined in Eq. ( I33p . the phenomenological model may 
be used to predict other experimentally-measurable quantities. In this section, we use 
the phenomenological model to predict the fraction of radiative events that satisfy 
requirements on the photon's energy and angular distance from the charged lepton. 
Radiative fractions for K® 3 decays were originally proposed by Ref. [23] and later 
refined in Refs. SI] . 

The radiative fractions are defined by the equation 

K° Km (E?», C n ) = T(^° - 7r-*Vy; E, > E™, h > 0^)/T K « , (35) 



f:i 



where E™ m is the minimum energy of the radiated photon, and 9J^ n is the minimum 
angular separation between the photon and the charged lepton, and T^o is the 
K® 3 decay rate including radiative corrections. Both E™ in and Qf" 1 are quantities 
defined in the kaon center of mass. 

In Table El we present the calculated radiative fractions. For the K Q e3 and the 
K® 3 decay modes we consider two values of E™ m : 10 MeV and 30 MeV. We do not 
make any requirement on the angular separation between the photon and the muon; 
however, for the K® 3 decay mode, we also consider 6^ n = 20°. The uncertainties 
in Table [3] are composed of two components: the uncertainty on the next-to-leading 
order correction and the uncertainty on the K® 3 decay rates. For the K® 3 and the 
K° 3 decay modes, we take the uncertainty on the next-order corrections to be 1.3% 
and 1.9% of T(K° -► n-e + v T , E 1 > Ef a , 6 h > 6^ in ), respectively. We estimate the 
uncertainties on the K® 3 decay rate by using the 0.3% uncertainty on the radiative 
correction parameters, S LD . 

6 Conclusions 

We have studied radiative corrections to the K® 3 and the K® 3 decay modes using 
a phenomenological model for their interactions. In our calculation, we include the 
contribution to the radiative correction from outside the three-body K® 3 Dalitz region. 
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Table 3: Radiative fraction, lZK i3 (E™ m , 6*™ n ), as defined in Eq. (1531) . Radiative frac- 
tions are given in percent (%) and they are presented for both the K® 3 and K^ 3 decay 
modes. 



n e3 (%) n° Kia (%) 

ffmin J gmin qo 20° 0° 

10 MeV 4.93 ±0.06 1.89 ± 0.02 0.564 ± 0.01 

30 MeV 2.36 ±0.03 0.956 ± 0.01 0.214 ± 0.004 



Using the cutoff scale A = m p , we find that 8 e LD = (1.3 ± 0.3)% and 5^ D = (1.9 ± 
0.3)%. The uncertainty in the long-distance radiative corrections is dominated by 
their sensitivity to the cutoff scale A and virtual calculation method. Dependence of 
the radiative corrections on the variation of the hadronic K-ir form factors is found 
to be small (0.04%). The calculated value of the 5 e LD radiative correction parameter 
agrees with a recent calculation [8] using ChPT. 

Combining the long-distance contribution from the phenomenological model with 
the short-distance contribution from Ref. [17], the total radiative corrections to the 
K® 3 and the K® 3 modes are (3.5 ± 0.3)% and (4.1 ± 0.3)%. The uncertainty in 
the total radiative corrections are dominated by the uncertainty in the long-distance 
component. When measuring the V us CKM matrix element, the uncertainty from 
radiative corrections should be ~ 0.0003. This uncertainty corresponds to a small 
fraction of the "external uncertainty" (0.0021) quoted for the \V US \ measurement in 
Ref. H]. 

The program KL0R was written to numerically evaluate the radiative corrections 
and to generate Monte Carlo events. Monte Carlo events produced by KL0R may be 
used to understand experimental detection efficiencies. 
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A Loop Integrals 

Virtual corrections to the K® 3 decay mode require the evaluation of loop integrals. 
Below, we define the loop integrals for one-, two-, three-, and four- point functions [571 
158] |4"2] . The corresponding vector loop integrals and their expansion in terms of 
momentum four-vectors are also defined. In Sec. [5], we numerically evaluate the loop 
integrals using the cutoff regularization scheme. Therefore, implicit in our definition 
is a high-mass, cutoff scale, A. 
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The scalar one-point function is 



Mm 2 ) = - [d'k l (36) 

Z7T J (k z — to 2 ) 

where misa mass parameter. 

The scalar two-point function is 

Bo(pl,ml,mj) = ±J (fc 2_ m ^ )(( l+ pi) 2_ m ^, ( 37 ) 

where p± is a four-momentum and TOj (i = 1, 2) are mass parameters. The vector 
two-point function is 

B»(pi; ml, to 2 ) = — / -— 2 M 2Y = P^ B M\ m cn m i)> ( 38 ) 

where -E> M (pi; TOo,to 2 ) is proportional to pi ijU . 
The scalar three-point function is 

2 222 2\ ± ' "^ 



^.^ftjmo,^^) = -J {k ,_ mlm + Vi? _ mlm + V2? _ ml) (39) 

where Sy = (pi — Pj) 2 is the squared difference of two four- momenta. The vector 
three-point function is 

rl 4 kk 
. 2 2 2\ / " "''V 

C lx (p 1 ,p 2 ;m ,m 1 ,m ^ - ' 



i7r 7 (A; 2 — Too) ((A; + £>i) 2 — m\){(k + p 2 ) 2 — to 2 ,) 
= Pi, M Ci+p 2 , M C 2 , (40) 

where C* = Ci(pj, s 12 ,pj] ml, ml, to 2 ,). 
The scalar four-point function is 

D (pl, s 12 , s 23 , Pi, p|, sis; ^o, to 2 , to 2 ,, to 2 .) 

= ±Jd 4 k[(k 2 -m 2 )((k + p 1 ) 2 -m 2 )]' 1 (41) 

x [((A; + p 2 f - to 2 ) ((fc + p 3 ) 2 - m^r 1 
and the vector four-point function is 

= ±-fd*kk»[(k 2 -ml)((k+p i y-ml)}- 1 

x [((k + p 2 ) 2 - to 2 )((/c + p 3 ) 2 - mj)]- 1 
= Pi,pDi + p 2 ^D 2 + P3,^D 3 , 

where A = D i (p 2 1 ,s 1 2,s 2 3,pl,p 2 .,s 13 ;m 2 ,m 2 1 ,ml,ml). 
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(42) 



B Method II Expressions 

As in Method I, the matrix element for the virtual corrections using Method II may 
be written in the following form: 

M v p l £ = -^V2G F V us f + (0) u(p v )P R [A pole i> K - B pole fa } v(p t ), (43) 

where A po i e and B voie are composed of loop integrals. In this section, we present the 
expressions for these parameters. 

We break-up the expression for A po i e and for B po i e into five components, 

5 5 

A po ie = ^2 Ai , B vo[e = s ^ j B i , (44) 

i=i j=i 

where the '1' component is from inter-particle photon exchange, the '2' component 
is from the 7r and £ wavefunction renormalizations, the '3' component is from photon 
exchange between the effective vertex and the charged lepton, the '4' component is 
from photon exchange between the effective vertex and the pion, and the '5' compo- 
nent comes from emission and re-absorption of a photon by the effective vertex. Since 
single photon radiation from the effective vertex is small, we expect that A§ and B§ 
(double photon radiation) are negligible [43] . 

The expressions for the Ai and the Bi may be written, compactly, in terms of a 
handful of two-, three-, and four- point functions. We define B^ = B^m 2 ; to 2 , A 2 ) 
and B\ = Bi(mj; mj, A 2 ) where i = 0,1. We also define the two-point functions that 
depend on the pole masses m + and m , 

B^ A) = B (t;\*,ml) 

B^ B) = B (t;ml,ml), (45) 

B^' C) = B (m 2 K ;ml,ml), 

where a = +, 0. We define the following four-point functions: 



C (a-,B) _ c i (s 7Te ,0,m 2 K ;ml,m 2 e ,ml) 
Cf> C) = aKV,0;m 2 ,A 2 ,m 2 ) 
C^ D) = a(m 2 ,*,m^;m 2 ,A 2 ,m 2 ), 



(46) 



where s^ = (Pi — Pj) 2 , i = 0, 1, 2 and a = +, 0. We define the following four-point 
functions: 

D\ o) = Di(ml, s K i, 0, t, mj, m 2 K ; A 2 , m\, mj, m 2 ), (47) 
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where Sy = (p, — Pj) 2 , i = 0, 1, 2, 3 and cr = +, 0. 
Using the above definitions, A\ and B\ are 

Ai = -2m 2 h Cj +,A) - 2AC B - ACj 4 + Sm 2 + x^D ( +) - 2X 2 AD 

+4m 2 + m 2 n D{ +) + 4x* t ADi + Am 2 + (m] + 2x^)/4 +) - 2m 2 , AL> 2 
+4m 2 (m 2 K - ml)D { + ] - 2tAD 3 , 



Bi 



AC A + ACf + AC 2 A - 2 AC B + [4x vt - 2X 2 ]AD + [4x^ - 2m 2 ] AD 1 



»(+) 



2_2 n (+) 



+ 4m|s rf Z^ + [4x rf - 2m|] AL> 2 + 4m 2 + m 2 K Df> + [Ax nl - 2t] AD 3 , (48 



where 



AC? 

A A: 



m 2 + C^ X) + (m 2 K -ml] 



c 



(+,*) 



C 



(0,X) 



m\D[ +) + (m 



K 



ml 



D\ +) - A (0) 



(49) 



and X = A, B. The components A 2 and B 2 are expressed in terms of the wavefunction 
renormalization constants [see Eq. (|26|) ]. 



A 2 = f + (t)l%)(6Z ir + 6Z i ) , S 2 = / 2 (t) f (<^ ff + ^ 



(50) 



where t = (px — p-w) 2 - The Aj and Sj components arising from photon emission from 
the vertex and absorption by either the pion or the lepton (i = 3, 4) are given by the 
equations 



A 3 
B 3 



2 n (+,c) 



/+(*) imjC^-ttCl 



?(0,A) 



r(+,C) 



and 

A 4 

B, 



f + (t) -2B?'*> + -^ABf' - 8x K „ C™ + C\ 



(A) 



mz 



r(+,C) 



ti+,0) 



c (+,c) 



2 mi 



a 



(0,C) 



ml 2 AC, 



+ 2t[c{°' C) -m; 2 AC{ C) ]} 
+/ 2 (t) {2fi(°' A) - 4B® - 2B[ e) + 2m 2 ct C) ~ 2tC[°' C) } 



(51) 



/+(*) {4^ 



(•^) 



*(£>) 



-2B, 



(C) 





25, 



to 



2(m 2 Y -m 2 )-i + A 5 



771 



^ D) } 



ACf } 



+AC ( o> - AC\ L 
AC { D) - ACf } + f 2 {t) {2B^ A) + B™ - B^ 



,c) 



2(m 



K 



777, 



t + A 2 - mi 



a 



(+,D)\ 
J 



+ -Mt) {-m 2 + B^ + AB^ + \m 2 + B^ 



+- 



7775 

1 r 

2 

1 



(B) 



2(777^-777 2 )-t+A 2 



G 



(+,-D) 



2(777 



K 



777, 



£ + A' 



•AP) 



act+<Q 



2 nf+.O) 



-AS, 



2^(0, £)) 



m^' 



(52) 
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Figure 8: Comparison of the PHOTOS (dots) and the KL0R (line) Monte Carlo gener- 
ators, (a) and (b) plot the cosine of the angle between the charged lepton and the 
photon (cos6*£ 7 ) for the K® 3 and the K® 3 decay modes, respectively. In plots (a) and 
(b) radiated photons are required to have an energy greater than 10 MeV in the kaon 
center of mass. Plots (c) and (d) compare the log of the radiated photon energy 
(log 10 (E 7 )) , where E 1 is in GeV. In PHOTOS the IR cutoff is taken to be 1 keV and in 
KL0R the photon mass was taken to be 1 eV/c 2 . 



C Comparison of KL0R to PHOTOS 

One tool often used by experimentalists to model radiative effects is a "universal" 
Monte Carlo generator called PHOTOS |10j . For the K® 3 decay mode, PHOTOS estimates 
the size of the radiative corrections in the leading-logarithmic (collinear) approxima- 
tion. Since it treats radiation from the charged lepton and the pion independently, 
PHOTOS should succeed in reproducing the photon energy distribution but fail to re- 
produce distributions of the angle between the photon and a charged particle. To 
illustrate this behavior, we generate the log(£' 7 ) and the cos 6^ distributions using 
PHOTOS and KL0R [see Fig. [8]. In K® 3 decays, PHOTOS significantly underestimates 
radiation in the central and backward directions, while for K° 3 decays, photon radi- 
ation is suppressed in the central region. For both the K® 3 and K® 3 decay modes, 
PHOTOS and KL0R agree for lower energy photons. The discrepancy at higher energies 
is most likely a result of the "soft photon" assumption in PHOTOS. 
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Though PHOTOS is a powerful tool, precision measurements of K® 3 decays require a 
more accurate understanding of radiative effects. With KLOR, we hope to satisfy this 
requirement and be able to study the structure of the radiative corrections. It should 
be noted that the results of PHOTOS can be improved if the K® 3 inner-bremsstrahlung 
matrix element is incorporated into the PHOTOS algorithm. 
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